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Abstract 

We analyze multi-matrix chain models. They can be considered as multi-component 
Toda lattice hierarchies subject to suitable coupling conditions. The extension of such 
models to include extra discrete states requires a weak form of integrability. The discrete 
states of the g-matrix model are organized in representations of sl q . We solve exactly the 
Gaussian-type models, of which we compute several all-genus correlators. Among the latter 
models one can classify also the discretized c = 1 string theory, which we revisit using Toda 
lattice hierarchy methods. Finally we analyze the topological field theory content of the 
2g-matrix models: we define primary fields (which are oo 9 ), metrics and structure constants 
and prove that they satisfy the axioms of topological field theories. We outline a possible 
method to extract interesting topological field theories with a finite number of primaries. 
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1 Introduction 



In this paper we intend to analyze matrix models made of g Hermitean N x N matrices with 
bilinear couplings between different matrices. Unless otherwise specified, by this we mean 
an open chain of q matrices, each linearly interacting with the nearest neighbours. These 
models have been already introduced and partially analyzed in (for other approaches to 
multi-matrix models, see 0,||>@]>[|1]j HI' El)- The reasons to go beyond two-matrix models 
are diverse. The extended two-matrix model provides a useful representation of c = 1 string 
theory at the self-dual point, ||; in particular it naturally incorporates the so-called discrete 
states, which appear to be organized in s?2 multiplets. We find it natural to ask ourselves 
whether such a construction can be generalized. The answer is affirmative: in the extended 
g-matrix model we do find discrete states organized according to representations of sl q . 
More recently it has been shown, ||, that c = 1 string theory at the self-dual point, i.e. the 
two-matrix model, is a huge topological field theory in which we can distinguish primaries, 
puncture operators and descendants. As we shall see, this holds for 2q matrix model too, 
although with new features (for example, the number of primaries is oo n ). 

On a more speculative ground one may remark that two-matrix models lead via hamil- 
tonian reduction to reduced models characterized by classical hierarchies ||,[10] which can 
be interpreted in terms of topological field theories coupled to topological gravity; in turn 
the latter can be put in correspondence with string vacua. The two-matrix model analysis 
suggests that, if we want to reach more interesting string or W-string vacua, we have to 
shift to matrix models with several matrices. Although we do not go as far as proving this, 
nevertheless many elements we find seem to support such a conjecture. 

Finally, to end the list of the reasons of interest on a more formal ground, we recall that 
the integrable hierarchy characterizing two-matrix models is the discrete Toda hierarchy, 
while the discrete integrable hierarchy characterizing multi-matrix models is a generalization 
of the latter. As we have already pointed out, we can extend these models by introducing 
additional (extra) states and couplings. While the extended two-matrix model does not 
present any essentially new features, the g-matrix models with q > 2 do. In fact we can 
have in general only a weak form of integrability of the extra flows (as opposed to the strong 
integrability of the ordinary cases). This form of integrability is nevertheless sufficient for 
all our purposes. 

The paper is organized as follows. In section 2 we review, mostly from Q, the main results 
concerning multi-matrix models and derive the flows of the extended g-matrix models. In 
section 3 we solve the coupling condition of Gaussian g-matrix models. In section 4 we 
introduce the discrete states and discuss their group properties. We then compute several 
examples of correlators in 2g-matrix models. Section 5 is devoted to the topological field 
theory properties alluded to above. In section 6 we introduce a few simple examples of non- 
Gaussian matrix models. Finally section 7 is devoted to an analysis of the discretized ID 
string. The latter can in fact be envisaged as a chain matrix model with bilinear couplings. 
It is interesting to rederive known properties of c = 1 string in our formalism. Finally 
two Appendices are devoted to the W-constraints in g-matrix models and to an explicit 
computation, respectively. 
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2 Mult i— matrix models: general introduction 



We review here some general results concerning g-matrix models, O. The partition function 
of the (/-matrix model is given by 



Z N (t, c) = I dMxdMi . . . dM q e TrU (2.1 



where M\ , . . . , M q are Hermitian N x N matrices and 

q q-1 

M a M a+1 

a=l a=l 



with potentials 



Pa 

^« = V4X a = l,2...,q (2.2) 



r=l 



The paS are finite positive integers. 

We denote by J^ip ltP2 ....,p q the corresponding g-matrix model. It has become moreover 
customary to associate to the generic g-matrix model ( |2.1| ) the Dynkin diagram A q . Occa- 
sionally we will stick to this convention and speak abou nodes and links. 

We are interested in computing correlation functions (CF's) of the operators 

T a)k = trM k a 

and possibly of other composite operators (see below). For this reason we complete the 



above model by replacing (2.2) with the more general potentials 



V a = J2 t a ,rM r a , a = l,...q (2.3) 
r=l 

where t ar = t ar for r < p a . 

In other words we have embedded the original couplings i a>r into infinite sets of cou- 
plings. Therefore we have two types of couplings. The first type consists of those couplings 
(the barred ones) that define the model: they represent the true dynamical parameters of 
the theory; they are kept non-vanishing throughout the calculations. The second type en- 
compasses the remaining couplings, which are introduced only for computational purposes. 
In terms of ordinary field theory the former are analogous to the interaction couplings, while 
the latter correspond to external sources (coupled to composite operators). Any CF is ob- 
tained by differentiating lnZjy with respect to the couplings associated to the operators that 
appear in the correlator and then setting to zero (only) the external couplings. 

From now on we will not make any formal distinction between interacting and external 
couplings. Case by case we will specify which are the interaction couplings and which are 
the external ones. Finally, it is sometime convenient to consider iV on the same footing as 
the couplings and to set t Qj o = N. 

The most popular procedure to calculate the partition function consists of three steps 
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(i) . One integrates out the angular parts such that only the integrations over the eigenvalues 
are left, 

r q N 

Z N (t, c) = const / J] II dA Q , i A(Ai)e l/ A(A 9 ), (2.4) 

■> a =li=l 

where 

q N q-1 N 

U = ^ ^ V a (\ aji ) + ^2^2 C a ,a+lK,iK+l,i, (2.5) 
a=l i=l a=l i=l 

and A(Ai) and A(A g ) are Vandermonde determinants. 

(ii) . One introduces the orthogonal polynomials 

= A™ + lower powers, ^n(A g ) = A™ + lower powers 

which satisfy the orthogonality relations 

J d\i . . . dA 9 ^ n (Ai)e M r? m (A g ) = h n (t, c)5 nm (2.6) 

where 

q oo q—1 

(2.7) 

a=l r=l a=l 

(Hi). If one expands the Vandermonde determinants in terms of these orthogonal poly- 
nomials and using the orthogonality relation (|2.6| ), one can easily calculate the partition 
function 

N-l 

Z N (t,c) = const Nl Y[ hi (2.8) 

i=0 

Knowing the h(c, i)'s amounts to knowing the partition function, up to an ^-dependent 
constant. In turn the information concerning the h(c, t)'s can be encoded in suitable flow 
equations, subject to specific conditions, the coupling conditions. Before we come to that, 
however, we recall some necessary notations. 

For any matrix M, we define the conjugate A4 

M = H^MH, Hij = hAj, Mij = M jh M z (j) = Af^-j. 

As usual we introduce the natural gradation 

deg[Eij] = j - i, where (-Eij)fc,/ = h,kSj,l 

and, for any given matrix M, if all its non-zero elements have degrees in the interval [a, b], 
then we will simply write: M G [a,b\. Moreover M + will denote the upper triangular part 
of M (including the main diagonal), while M_ = M — M + . We will write 

N-l 

Tr(M) = J2 M a 
i=0 
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The latter operation will be referred to as taking the finite trace. 

Coupling conditions. 

First we introduce the Q-type matrices 



f q 

\\dX a ^ n {\i)e^\ a r] rn {\ q ) = Q nm (a)h m = Q mn (a)h n , a = l,...,q. (2.9) 

a=l 



Among them, Q(l), Q(q) are Jacobi matrices: their pure upper triangular part is 2+ = 
J2i Ei,i+i- We will need two P-type matrices, defined by 

f q d 

J II ^ a ( — en(Al))e^ m (A g ) = Pnm(l)h m (2.10) 



J dA 1 ...o(A^ n (A 1 )e^(— Vm{X q ))= P mn (q)h n (2.11) 



The matrices ([T^) we introduced above are not completely independent. More precisely 
all the Q(a) , s can be expressed in terms of only one of them and one matrix P. Expressing 
the trivial fact that the integral of the total derivative of the integrand in eq^^) with respect 
to A Q , 1 < q < q vanishes, we can easily derive the constraints or coupling conditions 

P(l) + V; + c 12 Q(2) = 0, (2.12a) 
c a ^ a Q(a - 1) + V' a + c a , 0+ iQ(a + 1) = 0, 2<a<g-l, (2.12b) 
c q _ hq Q(q - 1) + Vg + V{q) = 0. (2.12c) 

where we use the notation 

Pa 

V^ l = ^rt ol>r Q r ~ 1 (a), a = 1,2,..., 5 

r=l 

These conditions explicitly show that the Jacobi matrices depend on the choice of the 
potentials. In fact they completely determine the degrees of the matrices Q(a). A simple 
calculation shows that 



Q(a) e [-m a ,n a ], a = 1,2 



where 



and 



mi = ijp q - 1) . . . (P3 ~ 1)(P2 - 1) 

m a = (p q - l)(p q -i - 1) . . . (pa+i - 1), 2 < a < 

m„ = 1 



n\ = 1 

n a = (p a -i - 1) ... (p 2 - l)(pi - 1), 2 < a < 

n q = (p q -i - 1) ... (p 2 - l){pi - 1) 
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Throughout the paper we will refer to the following coordinatization of the Jacobi matrices 



mi 



Q(l)=/ + + £5>(*)£ M _ i , Q(q)=I + + J2J2 b tt)Ei,i-i (2.13) 
i 1=0 i 1=0 

and for the supplementary matrices 

Q(«)=E E T i a) <j) E i,i-h 2<a<q-l (2.14) 



i l=—n 



Flow equations 

The flow equations of the q-matrix model can be expressed by means of the following 
hierarchies of equations for the matrices Q(a). 

' Q{a) = [Q k + ((3), Q(pt)], l<P<a (2.15a) 



d 

dt/3,k 



Q(a) = [Q(a), Q k _(/3)}, a<[3<q (2.15b) 



These flows commute and define a multi-component Toda lattice hierarchy, [12|,[jy]. 
Reconstruction formulae. 

The coupling conditions and the flow equations allow us to calculate the matrix elements 
of Q{a). From the latter we can reconstruct the partition function as follows. We start from 
the following main formula 

' lnZ N {t, c) = Tr(Q r (a)) , 1 < a < q (2.16) 



dt a 

It is evident that, by means of the flow equations for Q(a), we can express all the derivatives 
of InZjv with respect to the couplings t a ^k (i.e. all the correlators) as finite traces of com- 
mutators of the Q(a)'s themselves. In other words, knowing the Q(a)'s, we can reconstruct 
the partition function (up to a constant depending only on N). In particular we can get 

g^bMte).^))^, !<<><<, (2.17) 

It was already noticed in (l| that this equation leads to the two-dimensional Toda lattice 
equation. 



2.1 Extended q— matrix models. 

It is important to be able to compute the correlators not only of the states considered 
above, but also of new states, the extra states. To this end we enlarge the (/-matrix model 
by introducing in the potential U new interaction terms, as follows. We change 

N 

i=l 6i,...,6 9 
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in ( |P| , |2.5|) , and, accordingly, 

A*-£ = £ ( 2 - 19 ) 
6i,...,6 8 

in (|2,7j). Henceforth a,,&i,Cj, ... will denote non-negative indices. 

We denote by X6 1) ...,6 g the state specified (classically) by J2iLi ■ ■ ■ It is clear that 
when bi = for all i ^ a, this state reduces to r Q) 6 a , while the corresponding coupling g 
boils down to t a ^ a ■ Moreover the previously introduced bilinear coupling c a ^ a+ \ is nothing 
but the above g when all the hi = except b a = b a +i = 1. 

All the couplings and states that do not appear in the original model ( |2.5| ) are called extra. 
Exactly as in the original g-matrix model, we can introduce orthogonal monic polynomials 
£ n (Ai) and rj m (X q ) and define the Q(a) matrices. This is parallel to what happens in the 
the extended two-matrix model, S. 

However, unlike the extended two-matrix model, in the extended g-matrix model, we 
cannot in general define flow equations in matrix form like eqs.( |2,15a] , |2.15b ). This is a 



remarkable difference between extended two- and g-matrix models (with g > 2), and, at 
first sight, seems to spoil integrability and any possibility of exact calculation of the CF's. 
Fortunately this is not the case. What one has to do is not to calculate the flows of the 
matrices Q(a), but the multiple derivatives w.r.t. the couplings of hiZjv, i.e. the multiple 
derivatives of h n , and express them in terms of matrices Q(a). One can verify that such 
'weak flows' commute, and thus integrability is preserved, although in a weak sense. 
The procedure is as follows. We first introduce two series of functions, [jl[, 

&Ht,\ a ) = J J] dXpUM)^- (2-20) 



and 



t\t,X a ) = J f[ dX p e^ Vm (X q ). (2.21) 



where 



(3=a+l 



a—1 oo a—1 

^ = EE t P,k^p + E c ft/3+iWi- 

/3=lk=l /3=1 
q oo <J— 1 

Va = E E t ^,k^P + E c /3,/3+l W+l- 
f3=a+l k=l (3=a 



Obviously we have 

ti\t,X 1 )=U*i), r l $(t,X q ) = Vm (X q ). 

but for other values of a one sees immediately that £( Q ) and r/W are not polynomials. But 
they satisfy the orthogonality relations 

dX a &\t,X a )e v ^r ) $(t,X a ) = 6 nm h n (t,e), l<a<q. (2.22) 
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Eq.(2.9) provides a definition of the Q(a) matrix in this basis 



j d\ a &\t,\ a )X a e v ^ V £\t,\ a ) = Q nm (a)h m (t,c), VI < a < q. (2.23) 

Therefore the spectral equations follow 

X a & ] = Q(a)C ia) , l<a<q. (2.24) 
A q t/( q ) = Q(a)rj( a \ l<a<q. (2.25) 

where £ a and r] a represent the infinite vectors with components £g , £f, . . . and r)Q , rjf , . . . , 
respectively. 

With these bases at hand one differentiates h n , i.e. ( |2.6[) for n = m, w.r.t the appropriate 
couplings and evaluate the results when the extra couplings vanish. The result contains 
derivatives of £ n and r] n w.r.t to the couplings, which in turn can be evaluated differentiating 
( [2.6D with n > m or n < m. Finally one can express the result in terms of elements of the 
matrices Q(a), by making use of the above defined bases and rf^. Inserting this into 
the expressions of the correlators, i.e. into the derivatives of InZ^ w.r.t. the appropriate 
couplings, one can express the latter in terms of finite traces of polynomials in the Q(a)'s. 

From now on, whenever it is not confusing, we use the simplified notation Q(a) = Q a . 

The 1-point CF is easily found to be given by 

< Xaw*, >= Tr(Q? ■ ■ ■ Q a q q ) (2-26) 
The derivation of the two point functions, by the above procedure, is as follows 

^ ^ a 2 ir^ n 

< Xa 1 ...a q Xb 1 ...b q >— 



n=0 og ai ...a q ogb 1 ...b q 



d 2 h ri 



- = fdX- ZnX$ 1 ...>% q Vn+ [ d\U\T +bl -K q+bq Vn + 

dga 1 ...a q dg bl ...b q J og bl ... bq ' J 



%i...b q 



Then, using 




n-l 



d 



n-l 



dgb!...b q m=0 
we obtain 



nin , 



l—0 q 



771=0 



QT---Q a q 9 )(Qi---Q b Q 



(2.27) 



Along the same lines, we get 



< Xa 1 ,...,a q Xb 1 ,...,b q Xc 1 ,...,c q >= Tr< Q l 



ai+6i+ci Qa q +b q +c q 



(2.28) 
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+ 



(QT-Q7 
QT-Q a q q 



Qi 



bi+ci Qb q +c q 



+ 



(q\ i+c \..q 



Q b 1 K..Ql")[Ql\..Q^) i + p. 



where p. (c.p.) means permutations (cyclic permutations) of the sets {01, a g }, {b±, b q } 
and {ci, c q }. The RHS's of both Q2.27D and p.28[ ) are symmetric under the exchange of 
the x operators. This property together with the fact that the RHS's can be written down 
in terms of the Q a 's, which are calculable, expresses what we call weak integrability. 



3 Coupling conditions in Gaussian models. 

There are several methods to solve matrix models. One is based on W -constraints (see 
Appendix A) and will be occasionally used also in this paper. The most powerful however 
consists of solving the coupling conditions to obtain explicit expressions of the Q(a) matrices 
and, then, inserting these into the expressions of the correlators. In this paper we will mostly 
consider ^-matrix models in which the interacting terms are at most quadratic. These 
models can be solved in general and the solutions can be expressed by means of very general 
formulas (see below), the reason being that the coupling conditions reduce to a system of 
linear equations in the Q(a)'s. We would like however to point out that much more general 
(than Gaussian) g-matrix models can in principle be exactly solved. The only trouble is 
that, when the potential terms are more than quadratic, the coupling conditions are non- 
linear equations in the Q(a)'s and we cannot find such compact formulas as in the Gaussian 
models but we have to proceed case by case. We will see later on examples of non-Gaussian 
models. For the time being let us concentrate on the Gaussian ones. They are sufficient to 
reveal the topological properties of the corresponding matrix models. 

Let us first introduce a more convenient notation for Gaussian models. The q-matrix 
models with quadratic potential have the partition function of the form: 

„ q g 9-1 

Z= I ]JdM a eMj2(^ M a + u a M a ) + J2 c a M a M a+1 ) (3.1) 

a=l a=l a=l 

We notice that the linear terms can be eliminated by suitable redefinitions of the matrices 
M a . However it is often useful to keep them distinct (for example to study the topological 
field theory properties). Therefore, whenever this does not complicate the formulas too 
much, we will keep the linear terms. 

We will solve first the coupling conditions of 3- and 4-matrix models, both for pedagogi- 
cal reasons and in order to have explicit formulas of the simplest cases, and then the general 
case. 



3.1 The 3-matrix model 

The coupling conditions are : 

Pi +2tiQi + ui + ciQ 2 = 
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2* 2 Q 2 + u 2 + ciQi + c 2 Q 3 = (3.2) 
P3 +2* 3 Q 3 + n 3 + c 2 Q 2 = 

Eliminating the matrix Q we get the following two-matrix model type coupling conditions: 

U2C1. C\C 2 



P 



P 3 



+2(ti 
+2(* 3 



r- 2 

4* 2 

r- 2 



)Qi + («i 

)Q3 + («3 



2* 2 ' 
u 2 c 2 . 



2t 2 

ClC 2 



Q3 

Qi 











4t 2 ' v 2t 2 ' 2* 2 

Solving the system we obtain the following form of the Q matrices (with reference to the 
coordinatization ( |2.13 , 2~l~i| ) 

, 2t 2 

b = "-g 



T, 



(2) 



2t 2 , , 



CiC 2 Ui - 2c 2 t\U 2 + (4*1*2 - cf )«3 

2tit 3 u 2 + c 2 *iu 3 ) 



«0 



§ ((4*2*3 



cf )m - 2ci* 3 n 2 + cic 2 m 3 



and 



61 = -^(2*2(4*1*2-^ 



2? 



(2) 



P 

4ci* 2 t 3 n 



,(2) 



5 



2*in 
ci 



P 
P 



(2t 2 (4* 2 t 3 -cf 

2ciC 2 *2"- 



P 



where 



B = (4*2*3 - cl) (4*2*1 - c\) - (cic 2 ) 



while all the other coordinates vanish. 

So far we have used a basis £ n corresponding to the first matrix or first node and to 
a basis r] n corresponding to the third matrix or node. One may wonder what happens if 
one switches from the node 1,3 to the nodes 1,2. The coupling constraints are of course 
modified: V3 disappears from the third eq. (|3.3|) and in the second eq.( |3.3D there appears 
V 2 .Eliminating now the matrix Q 3 we obtain the 2-matrix model coupling conditions: 



T~2 



+2*iQi + Ul 

+2(* 2 



ciQ 2 = 



^)Q2 + («2 
4*3 



u 3 c 2 , 

2*3 ' 



ciQi = 



Calculating the form of the Q matrices we get the same result as above. Hence, changing 
the basis does not modify the model. This is the simplest example of a base independence 
property which must of course hold for all multi-matrix models with generic potentials. 



3.1.1 The 4-matrix model 

For the 4-matrix model the coupling conditions are: 

Pi + 2*iQi + ux + ciQ 2 = 

2* 2 Q 2 + u 2 + ciQi + C2Q3 = 

2*3<?3 +U3 + C4Q4 + C2Q2 = 

Vl + 2* 4 Q 4 + u 4 + C3Q2 = 



(3.3) 
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Eliminating the Q 2 and Q3 matrices we get the following constraints: 

4t 2 i3 - C 2 4t 2 *3 - C 2 4*2*3 - C2 

tT ,o/t c|* 2 , , 2* 2 n 3 ~ u 2 c 2 . C1C2C3 

P 4 +2(*4 - TT - : 2JQ4 + «4 - C 3 — — 2— - — 2Q1 = 

4*2*3 - C2 4*2*3 - C2 4*2*3 - «2 

For Q matrices we obtain the following form: 

_ s l ™(2) _ S 2 



-2) 

~ "A* T ° - A(* 2 *3 - ^ ' ) 



where 



A = 4*1*2*3*4 - 4c?* 3 *4 - 4c^*i*4 - 4c|*i* 2 + (cic 3 ) 2 (3.5) 



while we do not write down the explicit expressions for the s Q 's; they are linear functions of 
u a , therefore when u a = the Q(a) matrices are traceless. The other coordinates are: 



(2) 2*jn ( 3 ) _ 4*!* 2 - cf _ cic 2 c 3 n 

i_j — , i_i —n , /14 — — 

ci cic 2 A 



and 



2n(cf* 2 + cf* 4 - 4*2*3*4) 



A 

r.2 



^(2) cire(4* 3 *4 - c%) (3) _ 2cic 2 * 4 ra 

1 " A ' A 

, 2n(cf* 3 + ci*i -4* 2 * 3 ti) 

61 = A • 



Let us consider the two previous models at the cosmological point, i.e. when all the 
couplings are set to zero except the bilinear ones (the c a 's). The reason of the name is 
that in such a case the CF's essentially depend only on N, which is interpreted as the 
renormalized cosmological constant (see section 5). 

We see immediately that, while such a point is well-defined for the 4-matrix model, it 
is singular for the 3-matrix model (in fact A 7^ but B = 0). These two models reveal the 
difference between odd and even g-matrix models. The cosmological point is well-defined 
only for even g-matrix models. 

3.2 Gaussian q— matrix models 

Let us now concentrate on the most general case ( |3.1|) . In particular \x takes the form 

q q q-1 

H = Ai(Ai, ... ,X q ) = U a^a + ^2 ta ^ + X! c aKK+l (3.6) 
a=l a=l a=l 
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The coupling conditions are 



P(l) + ui + 2tiQ(l) + ciQ(2) = 

u a + 2i a Q(a) + CaQ(i + 1) + c Q _iQ(a - 1) = 0, a = 2, 
P{q) + u q + 2t g Q{q) + c 9 _iQ(g - 1) = 



(3.7a) 
(3.7b) 
(3.7c) 



These coupling conditions imply that Q(a) has only three non-vanishing diagonal lines, the 
main diagonal and the two adjacent lines. Now let us simplify the coordinatization of such 
matrix as follows 



Q(a) = e+(a) + e (a) + e_(a) 



(3.8) 



where 



e -(") =^29a(n)E„ jn -i, e (a) = ^s Q (n)P nin , e+(a) = ^ h a (n)E Hin+1 



with the understanding that h\(n) = 1 and <? ? (n) = R(n). In terms of these coordinates the 
above coupling equations take the form of the following linear system 



2ii + ci/i 2 (n) = 

2iisi(n) + cis 2 (n) + u\ = 

n + 2t 1 g 1 (n) + cig 2 (n) = 

2t a h a (n) + c a h a+ i(n) + c Q _i/t Q _i(n) = 0, a = 2, . . . 
2t Q s Q (n) + Co,s Q+ i(n) + c Q _is a _i(n) + u a = 0, a = 2, . 



2t a g a (n) + c Q 5 Q+ i(n) + c Q _i5Q,_i(n) = 0, 
n + 1 

A( TC + 1 ) + 2 *A( n ) + c «-iVi( n ) = 
2t q s q (n) + c g _is 9 _i(n) = 
2t q R(n) + c g _i5 g _i(n) = 



a 



2,. 



(3.9a) 



(3.9b) 



(3.9c) 



The solution of this system is expressed in terms of the matrices X a and Y a , defined as 
follows 



(3.10) 



and 





/ 2h 


Cl 


... 





\ 




Cl 


2t 2 


c 2 . . . 








X a - 







c 2 


2t 3 ... 
















... 


2i Q _i 


Cq-1 




V o 





... 


Cq-1 
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(3.11) 
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Of course Y\ = X q . One finds 



h a (n) = (-l) a (ciC 2 . . . Ca-l) ^etXa-i 

-1 



R(n) = (-l)'ncica • • • c q -i [BetXqj (3.12) 
SaW = (-1) ncic 2 ■ ■ ■ c a _i 

Moreover, if we denote by S and U the vectors (si, S2> . . . , s q ) 1 and («i, . . . , it ? )*, respectively, 
we have 

S = -X^U (3.13) 

As we have already remarked we can always without loss of generality suppress the linear 
terms in u a by constant shifts of M a . In such a case S = 0. 

It is now easy to see that, at the cosmological point (t a = u a = 0), the solution ( p. 12 ) 



is well defined when q is even, while it is singular when q is odd - in the latter case, for 
example, DetX q = 0. 



In the last part of this section we would like to dispel a seemingly obvious objection to 
the very content of this paper. Take the generic quadratic model of q matrices with nearest 
neighbour interactions 

q 5-1 q 

U=J2 ta M a + c a , a+ iM a M a+1 = M a A aP M p . (3.14) 

a=l a=l a,f3 

The q x q matrix A is symmetric, and, for the theory of central quadrics, it can be brought 
to a canonical diagonal form with all ones or minus ones on the diagonal. The signature of 
A is of course a characteristic of the potential. 

Let us see the consequences of this simple remark as far as the matrix model is concerned. 
The diagonalization of A can be achieved by integrating in the path integral over suitable 
linear combinations of the matrices M a , instead of integrating simply over the M a 's. Of 
course this gives rise to a Jacobian factor, which is however one if one uses only shifts of the 
Ma's. In this way one brings A to the diagonal form 

A = Diag(/ 1 ,...,/ g ) (3.15) 

but does not rescale its elements to unity. However this form is sufficient for our subsequent 
discussion. The initial matrix model appears at this point to be equivalent to the decoupled 
model with potential 

U' = Y,faM 2 a . 

a 

with partition function Z = Const(N)(fif2-..f q )~ N2 ^ 2 ■ We remark however that this pro- 
cedure is of no help if one has to compute correlation functions of composite operators, in 
that it screws up the definition of the states and renders the calculation of the correlators 
practically impossible. The procedure followed in this paper, i.e. the use of the generalized 
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Toda lattice hierarchy, has precisely the virtue that it allows the calculation of the exact 
correlators of significant composite operators. 

Finally let us remark that we can easily generalize the results of this subsection to the 



cases when in the potential are present, beside the terms of ( 3.14 ), also interactions of the 
type c a> pD a Dp where D a = Diag M a and /3 ^ a — 1, a, a + 1. In such cases the method is 
the same as in the chain models, with the only difference that the matrices X a and Y a will 
have, at the position (a,/3), additional non-vanishing entries c Q/ g if the latter are present in 
the potential. 

4 Correlation function of discrete sl q states 

In the previous section we have shown how to solve the coupling conditions of a given model. 
In this section we show how to calculate various correlation functions of composite operators 
(or discrete states). To start with let us illustrate a basic property of the latter: in the q 
matrix model they are organized in finite dimensional representations of sl q . 

4.1 sl q symmetry of the discrete states. 

We have shown in section 2.1 that we can enlarge the g-matrix model by introducing in the 
potential terms of the form: 

Q 

9a u -a q If D * a , with D a = Diag(M a ) 
a=l 

We call discrete states the operators Xai,...a q coupled to g ai ,...a q - We introduce also Xo,...o = Q 
as the operator coupled to go,...o = N. Classically, Xa ly ...a q is represented by J2k=i • • • K^k- 
These states carry a built-in sL structure. To see this one has to consider the following 
generators 



1 N ( 

k=o \ 



d d \ 

Xi+1 ' k fT\ ' 1 < 1 <Q~ 1 



8 _ a 



Hi form the Cartan subalgebra of sl g , while Ef- and E^- are, respectively, the raising and 
lowering operators of the Lie algebra sl q , corresponding to the roots: 

in the standard notation. The action on the states is as follows: 
1 



HiXa\...a q — n ( a i a i+l)Xa 1 ...a q i ^i jXai,...,ai,...,aj,...,a q — X 



i,jXai,...,ai,...,a,j,...,a q Xai,...,a,i±l,...,aj^l,...,a q 



Therefore the set {Xai,...a q = Yli=i ^Tk • • • Ki \i J2i=i a i = n } form an (unnormalized) 



representation of this algebra of dimension 



,k ' ' ' q,ki 

n + q — 1 
n 
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Although everything we do here can be repeated for g-matrix model with q odd, we 
concentrate from now on on the far more interesting case of even q. The main reason for 
this is the well-definedness of the cosmological point when q is even. This will allow us 
to give an unambiguous topological field theory interpretation of the corresponding matrix 
models, while such an interpretation does not seem to be possible for odd q. Therefore, from 
now on, unless otherwise specified we consider 2g-matrix models. 



4.2 General properties of correlators 

The correlation functions of the extended multi-matrix model are in general defined by 

d 3 

< X„M „M • • • X„(«) „(«) >= a • • • a " iZN 

a i '- a 2g a i '- a 2g og m (i) dg („) ( „) 

"l '■■■ w 2q "l '■■■ w 2q 

Our purpose in this section is to calculate the correlation functions in two simple special 
cases: the pure chain models where we set g ai ,...a 2q = except for go...a a a a+1 ...o = c a and the 
quadratic models where we have also the following nonzero coupling constants go...a a =2...o = 
ton 9o...a a =i...o = u a . As a consequence the CF's will be functions of c a ,t a ,u a and N. The 
chain models were referred to above as the cosmological point of the relevant 2c/-mm, while 
the quadratic models can be considered as quadratic perturbations of the latter. This second 
terminology is related to the topological field theory interpretation of section 5. 

To see some general properties of the CF's, it is convenient to use the IF-constraints 
(see Appendix A). We write down the W constraints in terms of them and obtain a set of 
(overdetermined) algebraic equations which in general one can solve recursively. 

The CF's, in the chain models, have the following symmetry property: 

< X (1) (1) (1) (1) • • • X (") («) (») («) > 

a) ,...a ...a. ...a, a) '...a. ...a. ...ai. 

1 ' l j 2q 1 ' i j 2q 

= < X (1) (1) (1) (1) • • • X (") («) («) («) > 

/L a, ,...a v ' ...a K ' ...ai, a\ ' ...a- ' ...a. ' ...ax ' 

1 J I 2q 1 ' ] I 2q 

This is due to the symmetry of the W constraints and to the invariance of the chain models 
under the exchange i <-> j. 

In the chain models the CF's satisfy (charge conservation): 

E [(41i + • • • 41i) - (42 + • • • 4 n J)i < x a(1 ) a d) • • • x w a w >= o (4.1) 

To prove the last statement we rewrite the constraint as follows: 

X Ol5oi...02, < Xoi.-oa, > +2^(^ + 1) = 

ai>l,a2>0...a25>0 

E a a g ai ... a2q < X01...02, >= ° 5 2 < a < 2g - 1 



ai>,...a a >l...a 2( j>0 

E 

ai>0,a2 9 _i>0a2q>l 



E a^oi...a2, < X01...M, > +^(^ + 1) = 
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We differentiate these equations w.r.t. g (i) (i), . . . , g („) („) and set g ai ... a2a = except 

a l ...a 2q a 1 ...a 2q 

go...a a a a+1 o = c a . One gets 



n 



a \ < X n W „(i) • • • X„(«) „(«) > +ci < Xno...oX„ (i) ■ • • X„(«) „(«) >= 

a l — a 2q a l ■•• a 2 9 a l •■■ a 2, 3 a l •■■ a 2q 



n 



X] a « } < „m • • • X„w „(») > + < (co-ixo...i,aa=i...o + 

CaXo...oc,=i,i-o)x„(i) „(i)---X„(«) „(n)>=0, 2<a<2q-l 

a 1 —a 2q aj ...a 2g 

ra 

E(fe) 

°2<7 < ^("...af) • • • *a< n >...«W > +C2 <? < X0...011X (D JD • • • An) >= 

, _^ 1 2q 1 2q 1 2q 1 2q 

Subtracting the even equations from the odd ones we obtain the result. 

The last property partially reflects the sfoq structure of the discrete states as it means, 
at the cosmological point, the conservation of the eigenvalue of H = Hi + H3 H h H.*iq-\- 



In the remaining part of this section we are going to compute exact correlators, i.e. all- 
genus expressions, from which we can extract the genus by genus expansion. To obtain such 
an expansion it is of fundamental importance that we can assign to each coupling a degree, 
denoted [•], as follows 

[g ai ,...,a q } =y + J2 i N ] = y ( 4 - 2 ) 

1=1 

and for each quantity such as free energies, correlators, etc., we can define a genus expansion, 
each genus contribution having a definite degree, 

00 

F = \nZ N , F=J^ F h, [F h ]= 2(1 -h)y (4.3) 

h=0 



In eqs^d^O), y,yi are arbitrary nonvanishing real numbers. 



4.3 The 1-point correlation functions 

To find explicit expressions of the correlators it is more convenient to switch from the W— 
constraints to the method based on the solution of the coupling conditions. 

4.3.1 Pure chain models 

We specialize here to the case in which all the couplings, except the bilinear ones, are turned 
off: 

U,a = 0. (4.4) 

In this case the coupling conditions take the form: 

P 1+Cl Q 2 = 0, (4.5) 
c Q -iQ Q -i + c a Q Q+ i = (a = 2, ... ,2q - 1) 

V2q + C 2q Q2q-l = 
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This linear system is so simple that we do not need to rely on the formulas of the previous 
section. We note that Q\ has only the first upper diagonal, and P%, which represents a 
derivative, has only the first lower one 



Qi 



fo 1 



V 



\ 



Pi 



( o 
1 



V 



/ 



The first coupling condition (|4.5|) gives now Q 2 = — l/cie_. Using the second for a = 3 
one finds then Q4, and so on for all the even Q's up to Q 2(? . For the odd ones the procedure 
is the same starting from Q\. 



+1 



( _ 1)fc £l 

C2 



Q2k = (-1 



Cl 



C2k 
1 

( 

C 2 fc-1 



(4.6) 



Now we come to the correlation functions which are expressed in terms of the Q-matrices 
by means of the formula: 



(4.7) 



<Xa lt ...,a 2q >=Tt(Q^Q 2 ^ ...Q%>). 

Due to the particular form of the Q's it is immediate to verify the conservation law ( [4.1|) . 
In order to have nonvanishing trace the number of 7 + and 7_ must be the same, i.e. 



«1 + 03 + . . . + Cl2q-l —0,2—0,4 



02q = 



The result for the one point functions can be found by means of the identities [7+, e_] = 7o 
and 7 + e_ = n (7q being the identity matrix and n = diag(l, 2, 3, ...)): 



<Xa 1 ,...,a 2q > =J2n=x ( n + °1 ~ 0,2) {n + 0\ - a 2 + 1) . . . (n + 0\ - 1) 

(n + a\ - a 2 + a 3 - a 4 ) . . . (n + a\ - a 2 + a 3 - 1) 
(n + ai- a 2 + ■■■ + a 2q -i - a 2q ) . . . (n + a 2? - 1) 

where x = — min[ai — a 2 , ai — a 2 + a 3 — 04, . . . , a\ — a 2 + ... — a 2q = 0]. 



(4.8) 



4.3.2 Quadratic models 

We write down some special 1-point correlation functions in the quadratic models. The 
derivation can be found in Appendix B. 



<r > - Tr( Q r )= YY (~l) fc 2-M / N + I - k ), , a 
where h a ,g a are defined by eqs. (|3.12j ). 
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This is an all-genus expression. In order to extract the genus h contribution follow 
the above described recipe. In particular for the 1-point functions we have the following 
expansion 

oo 

< r a , k >=£< r Q , fc > h N 1+k - 2h 
h=0 

( N + r \ 

The only dependence on N comes from I I and we can expand it as follows 



N + r 
r + 1 



j2 N i+k-2h b2hik) (410) 

h=0 



Using the last relation we can extract the genus h contribution: 

^ Tr(nn y y {-l) k 2- k r\b 2h (l - fe + 1) t _ k+1 _ 2h 
<r a , r > h -Tr(Q a ) = Eg (r - 2 l)lk!(l - k)l " ' 

■{Kg a ) l s r - 21 

with: 

hh(n) = r i r 2 • • • r 2h, b (n) = 1 

l<ri...r 2 h<ra 

4.4 Two-point functions 
4.4.1 Pure chain models 

For the 2-point functions we have to use eq.( |2.27| ). As an example we calculate the correla- 
tion functions of the form: 

< XO...a a =r...OXO...a p =s...O > = < Ta.r^.s > 

The formula ( p. 27 ) becomes, in this case, 



< XO...a ct =r...OXO...a g =s...O > r 



TV([(<&)+,Q£]) (a<(3) 
Tr(-[(Q r Q )_,Q^]) (a>/3) 



We take first the case a=even, (Q a ~ e_) and /?=odd (Qg ~ i+) and the two point 
function is written as: 

< T a , r T/3 t s >= Tr([Q„, Qp]) (a > j3 and zero otherwise) (4-12) 

Here the number of Q a and Qp in each trace must be the same, to balance the 1+ and 
e_'s (remember that a and (5 have different parity) so that r and s are forced to be equal. 
The traces can be evaluated as above and one gets: 

<T a , r Tp,s> = {qaqp) r Tr[(e-) r , I+]5 r=s = 

/N-r N \ 

= (q a qp) r S r =s £-£ n(n + l)...(n + r-l) (4.13) 
\n=l n=lJ 
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where q, 



C2 C 4 



i 



Ci C3 c /3-2 
C2 C4 "'Cp-i ' 



Cl C3 "'Ca-1 ' VP 

When a is odd and [3 even the result is the same with r exchanged with s. When a and 



(5 are either both even or both odd the 2-point function identically vanishes. 
As an example: 

< n, r r 2 , s >= S rs (ElT - Eli) n(n + l)...(n + r - 1) 

< T5 >r T2, s >= because a = 2 < /3 = 5. 



4.4.2 Quadratic models 

The 2-point correlation functions is a very important quantity in matrix models because its 
singularity indicate the existence of critical points and its scaling near them evaluate the 
anomalous dimensions of corresponding operators. In our case 2-point correlation functions 
permits also direct the calculation of the metric for the associated topological model (when 
the puncture operator is Q = 

Using again the equation ( 2.27] ) we write down the two-point correlation functions: 



<r 1;r T atS >=Tr[(Q£) + ,(Q*)_] 
Using the form of the Q matrices (jO, 3.12| ) we get the result: 



(4.14) 



< T 1>r T a , s >= E E 

< 2/ < r 0<k<l 0<i<l-k 
< 21' < s < h> < I 1 l'-k' <j< 2{V - k') 

i + j = l + l' — k — k' 



r \ l s 
21 J \ 21' 

N-l 



n=0 



n + 21 - 2k - i 

i 



J 

n + j 
3 



n + 21' - 2k! - j 
j 



n + i 
i 



where /(a) is: 



f(a) = (giK/gJY^g^s^h 



i+k J' +1 r-2l h l' -ls-2l 



a = 1 ... i 



To calculate the 2-point correlation we needed the quantity Tr(/™e1 l I+e^). We use the 
fact that for (n > m): 



Fle r 



k=0 



k T n-m-k A{n,m) 



A 



(n,m) 



n\m\ 



k\(n — m + k)\(m — k)\ 



Using this sum we get : 



n\m\(q + k)\ 



J k\{n — m + k)\{m — k)\ 



N 
+ k+l 



J n+p,m+q 



(4.15) 



(4.16) 



The evaluation of higher genus contribution follows the same method we have followed 
at the calculation of 1-point correlation functions. The only dependence of comes from 



sr^N-l P+s 



p + n 



and we are looking for the contribution of order N 



l+k-2h 
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We define the function B(r, s\n,m) as the coefficient of the genus h expansion: 

£j( P V ) ( ) =Nm+8 ~ 2hB h(r,s\n,m), r + s = n + m (4.17) 

The explicit expression is: 



if s + s'-l 
2 2h-l 



(r - s y h - l + 



(r _ s )i+2fe-i / s + s ; -Z + l\ 
+ I 1 - 2h - I 



B 2t s + s'-l \ s + s'-l-2t + l \ )2 h-l-2t+l 
,2<2 t <W-Z -2F\ 2t _ 1 2h-l-2t + l [r S) 



(r + r' = s + s'); where B^t are Bernoulli numbers and 7 are: 

I 

Ji(s,r,s') = H h---h3k+i---ji 

k=0 0<i 1 ...i k <s-l n-m<j k+1 ...j l <n-l 

The genus h contribution is: 

< Ti,rT a , s >h= X! X] X! 

0<2l <r 0<k<l < % < I - k 

< 21' < s < k' < I 1 l'-k'<j< 2(1' - kf) 

i + j = l + l<-k-k' 

rlSl \ L ) 1 tt„\ tktW -k-k '-2h 



-f(a)N l 



(r - 2/)!(s - 2l')\(2l -2k- i)\{2V - 2k' - j)\k\k'\i\jV 
(B h (2(i + k - I) + j,j\2l -2k-i,i)- B h (2(j + k' - I') + i, i\2l' - 2k' - 

5 Topological field theory properties of 2q— matrix 
models 



We study in this section the content of 2q-matrix models in terms topological field theo- 
ries. The motivation is offered by the example of 2-matrix model, which can be interpreted 
as a topological field theory with an infinite number of primary fields, @. We want to sec 
whether a similar conclusion can be drawn also for multi-matrix models. The easiest way to 
identify a possible topological field theory (TFT) content is to go to the cosmological point. 
We have seen previously that such a point is not well defined for odd q multi-matrix models. 
Consequently, in this section, we concentrate on even q multi-matrix models. To be definite 
we start with the 4-matrix model. We recall that the cosmological point is identified by 
setting all the couplings to zero except the bilinear ones, c a ^ a+ \, with reference to eq.fl2.ip. 
To simplify things further we set from now on 

C a ,a+l = ( — 1)° 
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without loss of generality (one can obtain the same results by suitable rescaling the couplings 
of the discrete states). Finally we replace N by a continuous variable x (i.e. we pass to a 
continuous formalism by suitably rescaling all the quantities and taking N —* oo; x is the 
renormalized quantity that replaces N). 

After these preliminaries let us concentrate on the 4-matrix models. Among the discrete 
states, our candidates for primary states are {ip a ,b,Q,^c,d}, where 

"0(1,6 = Xo,0,6,0, ^c,d = X0,c,0,d 

The relevant genus correlators to study the TFT properties can be computed from ( p. 27 ) 
and fl2~2§l ) 

< VWlVWa^d > = ( («1 + &l)( a 2 + h){c + d) - c(ai& 2 + Ml + &1&2) 



b±b 2 cd 



X C+d 1 ^ai+a2+6i+62,c+(i (5.1a) 



c + d-l. 

< ipa,bUc lt diUc2,d2 > = ((a + 6)(ci + di)(c 2 + d 2 ) - b(cid 2 + c 2 d\ + cic 2 ) 



abc\c 2 
a + b-1 



_a+6— 1 



"a+6,ci +C2+rfi +cfe 

(5.1b) 



and 

< Qi>a,bU c ,d >= (oc + ad + M)a; a+6_1 ^a+6, c +d (5.2) 

We will also need < QQQ >= x , which follows from the fact that the correlators involving 
only Q are the same as in the 2-matrix model, see [||. 

Now, at any level r = a + b let us select an arbitrary state among the ifi a b's and call 
it ip r , r > 0. Let us call C the collection of such choices for any r. Moreover, let us set 
u s = uq s . Then the states {ifj r ,Q,uj s } constitute the set of primary states of a TFT with 
puncture operator either Q or rpi or u±. This can be seen as follows. The non-vanishing 
structure constants are 

C n ,r 2 ,s =< ipnAzUs >= rir 2 sx s ~ l 5 S)ri+r2 

C r ,si,s 2 =< tprU sl U S2 >= rSiS 2 X r ~ 1 5 ryS1+S2 

Co, r ,s =< QA-u s >= rsx r " 1 5 r)S) C ,o,o =< QQQ >= x^ 1 

together with the ones obtained from these by permutation of the indices. Now, if the 
puncture operator is Q, the metric is 

rj r ,s = r]s,r =< Q*l>rU s >= rsx r ~ l 5 r , s , i]o,o = x~ 1 , (5.3) 

If the puncture operator is ipx, the metric is 

Vr,s = Vs,r =< ipiipr^s >= rsx r 5 s ^ r+ i, 7y 0) i = m,o = 1 ( 5 - 4 ) 

The case when the puncture operator is uj\ is exactly specular to the latter. These three cases, 
with exactly the same formulas for structure constants and metric, were met in ||, where it 
was proven that the inverse metric exists and the associativity conditions are satisfied. The 
TFT obtained with a definite choice C will be denoted Tq. If necessary one can specify the 
symbol of the relevant puncture operator. 
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Similarly, among the states oj C) d, c + d = s let us choose an arbitrary one and let us call 
it u 8 , s > 0. Let us call C such a choice for any level s. Moreover, let us set ip r fi = ip r , 
r > 0. Once again the states {tf) r ,Q,u) s } constitute the primary states of a TFT with 
puncture operator either Q or ipi or lo\. We do not need to explicitly prove this since the 
formulas for the structure constants and the metrics are the same as the previous ones with 
the substitutions ifj r — * tp r and uj s — * Q s . The TFT obtained with a definite choice C will be 
denoted Tg. 

We can think of Tq and Tg as unperturbed TFT's to which we couple topological gravity. 
Therefore we are going to have puncture equations and recursion relations. The latter are 
the same as in Q and will not be repeated here. The former can be derived from the W 
constraints. Instead of writing the most general formula, we write down the simplest one 
for the puncture operator tpi 

Xoi,a2,03,04 ^ — a 2 < Xoi,02— 1,03,04 

> +a 4 < Xai 

from which one can infer the action of the puncture: tpi lowers the even indices by 1. 
Therefore, when tpi is the puncture operator, the descendants of ip a ^ are going to be Xa,n,b,m 
for positive n and m, while any uj c ^ may be simultaneously primary and descendant, or an 
isolated primary. 

We notice that the situation here is an interesting generalization of the situation in 2- 
matrix model, pfl, where we have an infinity of primary states denoted {T n ,Q,T_ m }, with 
nonnegative integer n and m, where T n , T_ m are the discrete tachyonic states. Here we have 
oo 2 primary states, which depend on two integral indices and could be referred to as colored 
tachyons. In 2-matrix model, via reduction, one obtains an infinite set of TFT models (the 
n-KdV models) coupled to topological gravity, whose primary and descendants are to be 



found among the T n 's (or, symmetrically, among the Z!_ n 's), [11]. Similarly here we expect 
that, via reduction (see next section), the set ip a ^ with a and b positive, may support a 
series of matter TFT's coupled to topological gravity (i.e. primaries and descendants). Due 
to its characteristics - triangular structure of the primaries and relation with the product 
of two n-KdV models - possible candidates (certainly not the only ones) are, for example, 
the W3 topological minimal models coupled to topological gravity, ][D| , fi~9|l . 

In general, if we pass to 2g-matrix models, the set of primaries will be represented by 
the states {Xoi,o,a 3 A-,a2,-i,o}, by Q and by {xo,a 2 ,o,...,o,a 2q }- the primaries are oo q . This 
q should be related to the target space dimension in a string interpretation. In analogy 
with our previous discussion we are lead to speculate that one of the two sets above can 
accomodate the states of the W q+ \ minimal models coupled to topological gravity or anlogous 
TFT's. 



6 Non— Gaussian matrix models. 

The Gaussian version of g-matrix models is sufficient to study many properties, in particular 
it is enough to identify the TFT character of these models. From this point of view adding 
new interaction terms amounts to switching on new perturbations, which is not a very 
interesting complication in itself. However, if we come to reductions, i.e. if we try to extract 
TFT's with a finite number of primaries from the infinite TFT's that characterize the multi- 
matrix models, we are obliged to introduce non-Gaussian interaction terms. In this paper 
we limit ourselves to an example: our purpose is to show both the complexity inherent in 
non-Gaussian perturbations and a possible way to circumvent it. 
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The example consists of switching on a cubic potential in the 4-matrix model. More 
precisely we study the model ^3,2,2,2- The coupling conditions are the same as ( |3.4| ), except 
that in the first equation a term 3v\Q 2 must be added, where v\ = ii 3 is the coupling of 
the cubic term in the potential. Consequently the coupling conditions become a non-linear 
system of equations for the QaS. Eliminating Q2 and Q3 we obtain: 

d , nn c i ts \r> , 1 1 2t 3 n 2 -U3C2, C1C2C3 2 

4*2*3 - Q 4t 2 *3 - C2 4*2*3 - Q 

^ + " 4^3^! ^ + ^ " C3 4*2*3 -C| } " m^ 2 Ql = (6 - 5) 

We can therefore write this system in a simplified form as follows: 

Pi + 3viQ\ + 2*iQi + ui + cQ 4 + = 

V\ + 2* 4 Q 4 + u 4 - cQi = (6.6) 
These are formally the coupling conditions of the model ^3,2 with Q2 replaced by Q 4 , the 



couplings being suitably renormalized, |11]. 



With reference to the coordinatization ( 2.13| ) the equations (|6.6| ) in genus become 



2* 4 u 4 + ca (n) 

oi(n) = —R(n), bo(n) 



c 2* 4 
W = M") = -R(n) 2 



and the recursion relations 



2a (n) = -p- + -J-(c + ^1 (6.7) 



3ui 6* 4 ui ^ i?(n) 

2c* 1 c 3 \ c 2 u 4 cui 



2i?(n) = w a (nf + - r^-)ao(») - -g±- + 



2* 4 ^6* 4 t>i 12*4^1^ 12*4^1 6* 4 fi 

We have therefore to solve a cubic equation. Once we have done this all the unknowns can be 
determined and the explicit form of the matrices Q a can be calculated. The correlators (in 
genus 0) can be obtained as integrals of algebraic equations. Writing down such expressions 
is not very interesting. We can however ask ourselves whether in some region of the coupling 
space we can find some interesting solution. This is actually the case. 

Let us first simplify the formulas by imposing 6i>i = —1, t\ = CLj = 0. It can be shown, 
|p|], that this is no loss of generality. Then we impose the constraint 

a = (6.9) 

The above equations then imply 

ti , , , n n 2 



ai=ui, R(n) = --r, b = bi = 0, b 2 (n) = 



and, therefore, 



u\ = -^-n (6.10) 
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In general, when we impose constraints on the fields or on the coupling space, we are not 
allowed to use the flow equations of the integrable hierarchy (in this case the Toda lattice 
hierarchy), because such constraints might deform the dynamics in a non-integrable manner. 
However one can prove that, in the case of the constraint (|6.9| ) or, which is the same, ( |6.10|) , 
the constrained dynamics is still integrable and coincides with the Toda flows constrained 
by eq. (|6,9| ): the resulting hierarchy is the KdV hierarchy, [JO]]. In other words, if we look at 
our 4-mm in the submanifold of the parameter space specified by eq.( |6.10 ), the correlators 
are those of the KdV model. As is well known this model has only one primary field, ifti, 
and, in particular, < tpiipi >= a\ = u\. 

By considering the model M.p+2,2,2,2 one can find the p-th critical point of the KdV series. 
One can also identify higher KdV models. These are all (trivial) generalizations of the 2- 
mm. However, as we have remarked above, the structure of the q-matrix models allows for 
more complex and interesting reductions which will be the object of future research. 



7 The discretized ID string 



It is interesting to return to the problem of discretized ID string theory within the present 
formalism and see how we recover the results already obtained with other methods. 
Let us start by studying and solving the following Gaussian partition function: 



Tr {^Mi + c^MiMi+t 



2q 



2q-l 



2 z 



DMi exp 

The corresponding coupling conditions are {t = t/c): 

Pi + ctQx + cQ 2 = 0, 
tQi + Qi-i + Qi+i = 0, i = 2, 
T>2q + ctQ 2q + cQ 2q -i = 



(7.1) 



(7.2) 



,2q-l 



We express all Q a matrices a = 2, ... 2q — 1 in terms of Q\ and Q2 q - For this purpose we 
introduce the determinant ofuxn matrix: 



Dn 



t 1 ... 

1 t 1 

01 t '■■ 

. 1 '•• 

... '•. 

... 
The determinant D„ satisfies the recursion relation: 



1 

t 1 
1 t 



D n +i — tD n — D n _i 
D = l,D 1= t 



The expression of D n is : 



[n/2] I 
D n = E t n - 2k (-l) k 

k=0 \ 



n — k 
k 
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or: 



r n+l „n+l 
n - r 2 



where n, r 2 are the roots of the second order equation r 2 — tr + 1 = 0. 
Solving the system ( fflj| ) we get (Do = 1): 



Qi 



1) 



t+i 



D 2q - 2 



■{Q x D 2q ^ 1 - Q 2q D^ 2 ),i = 2, . . . 2 9 - 1 



For Qi,Q 2q we get the usual 2-matrix model with quadratic potential coupling conditions: 



Pi + ^rf 
The Qi,Q 2q matrices are: 



D 



2q-3 



D 2q - 2 



Qi + „ — Q 2q = o 
u 



D 



2q-3 



D 



2q-2 



2q-2 



Q2q + „ Ql = 

L>2q-2 



where: 



Qi = 1+ +gi(n)e-, 
Q2q = h 2q (n)I + + R 2q (n)e^ 



1+ — ^ En,n+1> e - — 2L n ^n,n-l, 



n=0 



n=0 



with 



and 



gi(n) = h 2q (n)R 2q (n) 



2 (ctP 2q - 2 - D 2 q--i)D 2 q- 2 ^ R ^ ^ ^ cD 2q _ 2 



B 



D 



B = c 2 (4i 2 - l)D 2 q _ 2 - 8ctD 2q - 2 D 2q - 3 + D 2 2q _ 3 
The form of Q a matrices is: 

Qa = h a (n)I + + g a (n)e-, a = 2...2q— 1 

with 

= — Z7=> — [cD 2q - a -! - (ctD 2q ^ 2 - D 2q - 3 )D a _ 2 ) 



9a(n) 



cB 
-l) a+1 
~B 



[{ctD 2q _ 2 - D 2q -z)D 2q - a _i - cD a _ 2 ] 



Using (2.26), ( 2.27| ), we can now calculate the 1-point correlation functions: 

(h a (n)g a (n)) r 



, T 2 ,_ Ti( 2 y (-l)* 2 -*(2r)l / iV 
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and 2-point correlation functions: 

<ri )2r r Q>2a >=Tr[(Qf) + ,(Q 2s )_] 
The 2-point correlation function is : 

< Ti,2rTa,2s >= E E L r3 (k,k' ,i,j\N)f(a) 

0<k<r 0<i<r-k 
< k' < s s — k' < j <2(s — k') 
i+j = r + s — k — k' 



with: 



20! (2^)! , ^k+k'o-k-k' 



L r , s (k,k',i,j\N) - {r _ 2l y {2l _ 2k _ i y kl{s _ 2m2s _ 2k ,_ j y kn ( ^ + 2 
A ! " / n + 21 - 2k \ ( n + 2i + 2j \ ( n + 21' - 2k' \ ( n + 2i + 2j 



E 



j 



n=0 

and where /(a) is: 

/(a) = (9ih a / 9a y +k (9a) r+s (h a y- r , a = 2 . . . 2q - 1 

We can obtain a more explicit formula for the determinant D n when 2 > t > — 2. In 
this case the roots of the equation r 2 — tr + 1 = are complex r± j2 = exp(±iu;/2) and the 
determinant is : 

sin(nu;/2) 
sin(io/2) 



with to = (l/2)arctan-\/(2/t) 2 — 1. This formula permits us to single out the dependence of 
the 1-point correlator < T a ^ r > on the parameter a. First we calculate: 

h a (n)g a (n) = (1/2)[(2A + A') + (A' - 2Acos(io{q - l)))cos(w((/ - a))} 

with 

= D 2q -3 ~ ctD 2q _ 2 = c 2 D 2 2q _ 2 - {ctD 2q -2 ~ D2g-3) 2 
B ' Be 

Hence we have the following behaviour: 

{K{n)g a {n)Y = £ d k e™ ka 

k=—r 

Now we pass from the discrete variable a = 1, . . . 2q to a continuous time t £ [0, T/2]. 
We introduce the puncture operator: 



C>( 2p ) = / dtTr(Q 2r (t))e ipt = K r V d k 5(p + uk) 



rt/2 

dtTr(Q zr (t))e tpc = 

where -fC r behaves like N r+1 . 
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In the pure chain models (no potentials V a ), Q(t) is proportional to either I + or J_. 
From this it follows that Tr(Q r (t)) is independent of t. Hence: 

(2p) ~ S(p) 

The conclusion is that for quadratic models we have, apart from the fundamental state with 
zero momentum, also other excited states (discrete states) with integer momenta p = nuj, n 
integer. 

We can study the 2-point correlation functions in the same framework. We take the 
particular case: 



< n,2rTa,2r >= E M n (g a ) n (h c 



\2r—n 
LV1 nKya) K n a) 

n=0 

with n = i + k and: 



M n = J2 L 2r , 2r (k,k',n-k,2r-k' -n\N)( 9l ) n 

0<k,k'<r 

and look at the dependence on the parameter a. Using the dependence on a for : 

(g a ) n (ha) 2r ~ n = £ ^ fc) exp(io;a) 

k=— r 

we can calculate the 2-point correlation function. Passing to the continuous time and using 

(k) (k) 

the symmetry of Ki- K_l we can write the expression of the 2-point correlation function, 
as: 

r 

< T\ : 2rT a ,2r > («) = E sinu;A;(a — a\) 

fc=0 

We can now evaluate 2-point correlators of puncture operators in the momentum space: 

r T/2 r K (k) 

G {2r \ P ) = / £e-**-*0 < n, 2r r a , 2r > (t) = E 2 f (7.3) 



We have intermediate states at all integer momenta p = kto, k integer. The pulsation 
to depends on the scaling we use when passing from discrete values of a = 1, . . . 2q to 
continuous time t E [0, T/2] . 

We can now apply all this to the c = 1 string theory model model, [16|[17|. The c = 1 
model with discrete time can be formulated as a multi-matrix model with the partition 
function: 



Z = J dMi exp 



Mi;«±^+ E x>-(M, 

\j=l i=l 



with a quartic potential V(M) = M 2 — gM 4 . However, only the contribution near saddle 
point V'(M C ) = 0, where the potential is quadratic in the fluctuation AM, is essential 
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The new partition function is (up to the constant exp(— N j3e/ {&g))\ 




It represents a string theory on circle with radius R ~ -. This is exactly our initial partition 
function ( |7.1| ) with the identifications i = 2(2e — 1/e), c = 1/e than t = 2(2e 2 — 1) and we 
have the following limiting cases for the determinant D n : 

D n ~ t n for e — > oo, t — > oo 

£>3n+fc ~ (-l) n (l + A:), fe = 0,-l,-2 for e -» 0, t -» 2(2e 2 - 1) 

From the results before we have found that our model describes particles with energy 
levels equal to (n + l/2)o;(e)//3. For small e, cj ~ e and changing e (lattice spacing) means a 
linear change of energy scale. In this limit the model describes a string with the discretized 
real line as target space. 

Because sin{ui/2) = 2e\/l — e 2 , for e < 1 the pulsation becomes uj complex, which is 
a sign of instability of the model. For e — > cxd the model decouples in g noninteracting 
1-matrix models. The instability is due to the liberation of the vortices and give rise to the 
Kosterlitz-Thouless transition for e near 1. 

Our method allows us to calculate (at least in principle) all n-point functions at any genus 
and what is more important it permits calculations in the vortex region where e < 1, t > 2. 



Appendix A. The W— constraints 

This Appendix is devoted to the derivation of the W-constraints in q-matrix models. 
From both the coupling equations ( |2.12cD and consistency conditions (2.15), we get the 
W-constraints in the form: Tr{Q n+r {a)d r Xa (*)) = where * are the relations (2.9) []. 
W-constraints have the form: 

W^ ] (a)Z N (t,c) = 0, r > 0,n > -r; a = l,...q. (7.5a) 

or 

(C^(a)-(-l) r T^(a))Z N (t,c) = 0. 

M 

involving the interaction operator Tn which depends only on all the couplings g ai ...a q , except 
50,...,0,aQ,,0,...,0 = ta,a a - 

For example T n ^ and T n 2 ^ are: 

d 

T n l \a) = a a g ai ... a - (7.5b) 

<jg a i...,a a +n,...a q 

T n 2] (a) = a a a' a g ai ... aq g a/i _ a ,— + 

°9a\ +a[ ...,a a +a' a +n,...a q +a' q 
9 

+ a a (a a - l)g ai . 



dgai...,a a +n,...a„ 



*For another approach see p(l 
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The operator Cn\l) has the same form as that of the two-matrix model: 

£M(1) = fdz: ^{d z + J) r+1 : z r+n (7.5c) 
where :: is the normal ordering and J(z) is the U(l) current: 

Pi OO o 

J(z) = J2 k^ k z k - 1 + Nz- 1 + £ (7.5d) 
k=i k=i dt ^ k 

The same expression holds for C n (q). 

The expression of £n\a), a = 2, ... q — 1 is different due to the absence of the P-matrix 



term: 



Cl\a) = Jdz: l -{d z + V>) r P a : z r+n (7.5e) 



with 

V^f^kt^z^+Nz- 1 , (7.5f) 



fc=i 

OO 



Pa = Nz- 1 + Y,Z- k - 1 1 £- 



The explicit expression of the first terms is: 

4V) = I>WoT + Wta,l<*n,-1 

fc=l 0*a,k+n fcljfc2 Ot a , k+n 

a, 1^71,-1 

As an example we write down the , and wj^ constraints for the three matrix 
model. 



WW: 



kt k ( T k-i) + N h + C12 (Ai) + Cis <<Tl) = 
^ <A fe _i) + iVui + ci 2 (n) + c 23 = (7.5g) 

E A;s fc (cr fe _i) + iVsi + c 23 (Ai) + ci 3 (ri) = 



iV(iV + 1) 



X! ( T k) + C 12 (Xlio) + Cl3 (XlOl/ 

Efcufc(A fc ) +ci 2 (xno) + c 23 (xon> = (7.5h) 
X! + C 13 (xioi) + c 23 (xon) 



S^fc<T5fc + i) + (JV + l)<Ti)+Ci2<X210)+Cl3<X20l) = 

ku k (Afe+i) + C12 (X120) + c 23 (X021) = (7.5i) 

E^K+l) + ( Ar + 1 )( CT l)+ C 13(Xl02)+C 23 (X02l) = 
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One easily sees from the second group of identities that the limit of pure chain models 
(cosmological point) does not exists for three-mm. The same thing holds for odd-q matrix 
models. However, writing down the W constraints for even-q matrix models, one can see 
that such a limit exists. This confirms the results obtained with other methods. 



Appendix B. Explicit derivation of lp correlators. 

In this Appendix we give the derivation of the 1-point functions promised in section 4. For 
this we need to know Q p where Q = 1+ + ae_ ;we express it in terms of the normal ordered 
quantities : Q l : (QQ =: QQ : +[QQ}) 

qp=Y,-- Q p ~ 2k ■ ([QQ]) fc 4 p) (7.5j) 

fc=0 

where the contractions are [QQ] = —clIq because [i+,e_] = Io ;the normal ordering means 
that we have expressions with I + on left side and e_ on the right side. 

The coefficient is the number of ways in which we can choose k pairs from p identical 
objects: 

a(p) _ 1 f p\ f P~ 2 \ fp~2k + 2\_ p\2- k 



k\ \ 2 J y 2 J y 2 J (p-2k)\k\ 
Using : 

. QP-2k ._ ^ I P~ 2k \ aii ,_2fc-i (e _ )j (? 5k) 



I 

i=0 \ 



We have the result: 

2p 2p~2k 



Q 2 v = J2 E \(-l) k a i+k A^ p) ll P ~- 2k -\e^y (7.51) 

k=o i=o V 1 ) 

We define the Q\ matrix: 

Qi = 1+ + a Io + ai(e-) = Q(ai) + a I 

Q i=P o ( 21 ) Q^ a o' 21 ( 7 - 5m ) 
Using (|7.5|) the expression of Q\ is: 

QI=EEE( 21 )( 2l ~ l 2k ) (-l) k 4 l) < +k ^- 2l ir 2k -^-Y (7-Sn) 

21=0 k=0 i=0 \ ) \ ) 

The same expressions are for : 

r 21 21 2k / \ / \ 

Q r a= EE E 07 (-^f/lW^)' (7.BO) 



21=0 k=0 i=0 



21 I \ i 

\i+k u 2l r-2l 



h{a) = {g a /K)^Ks r - M ,a = l... , 
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Because we have the summation: 

N-l 

I n 

k I \ k+1 



N=0 

the 1-point correlation function is: 

rr^r V- {-l) k 2~ k r\ (N + l-k\ n r _ 2l , 

<r r >=TVgi=^^ ? -^-4-— — (W-; a (7-5p) 
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